Sterile neutrinos: direct mixing effects versus induced 
mass matrix of active neutrinos 



Alexei Yu. Smirnov"''' ^ and Renata Zukanovich FunchaP'" ^ 

" ICTP, Strada Costiem 11, 340 U Trieste, Italy, 
^ Institute for Nuclear Research, Russian Academy of Sciences, Moscow, Russia 
'^Instituto de Fisica, Universidade de Sao Paulo, C. P. 66.318, 05315-970 Sao Paulo, 

Brazil 

Abstract 

Mixing of active neutrinos with sterile ones generate "induced" contributions to the 
mass matrix of active neutrinos ~ m5 sin^ 6as, where nig is the Majorana mass of the 
sterile neutrino and 9as is the active-sterile mixing angle. We study possible effects 
of the induced matrix which can modify substantially the implications of neutrino 
oscillation results. We have identified the regions of ms and sin^ 9as where the induced 
matrix (i) provides the dominant structures, (ii) gives the sub-dominant effects and 
(iii) where its effects can be neglected. The induced matrix can be responsible for 
peculiar properties of the lepton mixing and neutrino mass spectrum, in particular, 
it can generate the tri-bimaximal mixing. We update and discuss bounds on the 
induced masses from laboratory measurements, astrophysics and cosmology. We find 
that substantial impact of the induced matrix is possible if ms ~ (0.1 — 0.3) eV and 
s\T?9aS ~ 10"^ - 10^2 or > MeV and sin^ 6*^5 < 10"^. The bounds can be 
relaxed in cosmological scenarios with low reheating temperature, if sterile neutrinos 
decay sufficiently fast, or their masses change with time. 



1 Introduction 

There are two salient properties of neutrinos related to neutrality which distinguish them 

from other known fermions: 
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- the possibility to have a Majorana mass terms; 

- mix with new fermions which are singlets of the Standard Model (SM) symmetry group. 
So it would seem natural to explain unusual properties of neutrinos, such as smallness 

of masses and large mixing, using these two features. In fact, the see-saw mechanism [1] 
employs both. 

It may happen however that the see-saw mechanism is not enough to explain the pattern 
of the lepton mixing, especially if quark-lepton symmetry or unification are imposed. In this 
connection, we will concentrate on the second feature - the possibility of neutrinos to mix 
with singlets of the Standard Model, i.e., sterile neutrinos. In general, sterile neutrinos may 
originate from some other sectors of the theory, e.g. related to supersymmetry breaking or 
extra dimensions, and not coincide with the right handed components of neutrino fields. 

In all, there are two types of effects of the mixing: 

• Direct effects - when new states can be produced in various processes and participate 
in neutrino oscillations, etc.] 

• Indirect effects - via the modification of the mass matrix of light active neutrinos. 

The role and relevance of these effects is determined by masses and mixings of sterile 
neutrinos. In some ranges of parameters the cosmological and astrophysical consequences 
of mixing are more important and the influence on the mass matrix is negligible. In other 
regions vice versa: direct mixing effects are negligible - mixing becomes "invisible" but the 
modification of the mass matrix is substantial. 

Being light, sterile neutrinos can immediately take part in the phenomenology of neutrino 
oscillations changing the interpretation of experimental results [2, 3, 4]. Being heavy and 
weakly mixed, they do not show up in oscillations and other processes, however their "invis- 
ible" mixing can strongly modify the mass matrix of active neutrinos and therefore change 
implications of neutrino results for theory. In particular, in this way, the presence of sterile 
neutrinos can induce the large or maximal mixing of the active neutrino components [5, 6], 
or, for instance, produce deviation of the 1-2 mixing from maximal [7]. 
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In a range of masses ms < Q, where Q is the energy release in processes, sterile neu- 
trinos do not decouple since they still can be produced. They can decouple here in the 
sense that their direct dynamical effects are negligible due to the smallness of mixing. For 
bigger masses, ms > Q, sterile neutrinos decouple, as it happens in the standard see-saw 
mechanism. This decoupling generates an additional contribution to the mass matrix of 
active neutrinos and gives negligible deviation from universality. 

From observational point of view till now there is no clear evidences of existence of 
sterile neutrinos, though some hints exist. Those include, the LSND result [8] and its 
interpretation in terms of oscillations in the (3 -|- 1) or (3 -|- 2) neutrino mixing schemes; 
large scale structure formation in the Universe with the warm dark matter composed of the 
keV sterile neutrinos [9, 10]; high observed velocities of pulsars and their explanation as an 
asymmetric emission of the keV sterile neutrinos [11]; the early reionization of the universe 
due to the radiative S - decay [12]. 

In this paper we study in details the possible effects of sterile- active mixing on the 
mass matrix of active neutrinos. We obtain bounds on these induced masses from the 
direct mixing effects. We find the impact of the induced matrix may be considerable if 
ms - (0.1 - 0.3) eV and sin^ Oas 10"^ - lO'^ or ms > 300 MeV and sin^ Oas < 10"^. 

The paper is organized as follows. In sec. 2, we determine the mass matrix induced 
by mixing of active neutrinos with a sterile one. We study the properties of the induced 
mass matrix and the possibility to explain certain features of the neutrino mass spectrum 
and mixing pattern using this matrix. We find values of mixing (as a function of the 
mass of the sterile component) for which the effect of the sterile neutrino (i) explains the 
dominant structures of the mass matrix; (ii) produces the sub-dominant structures of the 
mass matrix; (iii) can be neglected, being of the order or below the la uncertainties of 
the present measurements. In sec. 3, we consider various bounds on masses and mixing 
of the sterile neutrino and consequently, on the induced matrix, in particular those from 
astrophysics and cosmology. We then, in sec. 4, confront these bounds with regions found 
in sec. 2 and discuss how they can be improved in the future. We also comment on new 
physics scenarios which allow to evade the bounds. Our conclusions are given in sec. 5. 
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2 Active-sterile mixing and induced mass matrix 



2.1 Induced mass matrix 



Let us consider three active Majorana neutrinos Va — (fg, J^t)'^ with mass matrix 



/ 



e 



^/ue '^^i^i '^jj.T 



\ 



TTlfQ TTlf^ TTlq-q- 

generated, e.g., by the see-saw mechanism [1]. We consider nia < 1 eV - below the present 
upper bound to avoid strong cancellations of different contributions. 

We assume that (i) the active neutrinos mix with a single (for simplicity) sterile neutrino, 
via the masses 

"^Is = i'^eS: m^s: ru^s) ] (2) 
(ii) S has a Majorana mass, ms, which is much larger than the mixing masses and ma- 
ms > ruas: rua. (3) 
So, in the basis {ua, S), the complete mass matrix has the form 



ma ruas 



(4) 

rns j 

Properties of S: masses, mixing, possible new symmetries, etc., arc determined by some 
new physics which, in general, differs from physics responsible for generation of ma- 

Under condition (3) the block diagonalization gives for the light neutrinos the mass 
matrix 

f5i ma + mi, (5) 

where 



mi = 



ms 



{mas) X {rnasf, 



(6) 



is the induced contribution to the neutrino mass matrix due to active - sterile mixing, or 
shortly, induced mass matrix. For the individual matrix element we have 

msmjs 



ms 



(7) 



Let us introduce the active - sterile mixing angles 

sine.s^'^. (8) 
ms 

Then the induced masses can be written as 

{mi)ij = - sin 9is sin 9jsms. (9) 

It is this combination of parameters which determines physical effects. For the flavor blind 
mixing we would have simply the product sin^ Osmg- 

In the case of a single sterile neutrino the induced contribution is the singular (rank-1) 
matrix. This feature substantially restricts possible effects of the induced matrix. In the 
case of two (several) sterile neutrinos, two (several) independent singular contributions to 
the induced matrix appear: 

i my 

That opens new possibilities in the description of neutrino mass matrices. Apparently, with 
three neutrinos any structure of the matrix can be reproduced. 

2.2 Neutrino mass matrix in flavor basis 

To evaluate the impact of the active-sterile mixing, we reconstruct the neutrino mass matrix 
from the data in the flavor basis in the context of three active neutrinos. The values of matrix 
elements in terms of the oscillation parameters are given by 

m^p = mi e-''^^ U^.U;, + K.U;, + m, e'^^^^ ui^u}^ (11) 

with a, P — e, /i, r, and Uai being the elements of the PMNS matrix. The matrix elements 
Uai are functions of the three mixing angles 9i2, Ou, 023 and the complex phase S given by 
the standard parametrization of the mixing matrix. 

We use the best fit values and 1 a intervals of the oscillation parameters from Ref. [13]: 

|Am^2l = 2.4 (l.OO to.ll) x 10"^ eV^ (12) 
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Am^i = 7.92 (1.00 ± 0.045) x 10"^ eV^ (13) 

sin2^23 = 0.44 (1.00 l°:?i), (14) 

sin^ ^12 = 0.314 (l.OOlH^s) , (15) 

and for sin^ ^13 we take 

sin^ ^13 = 0.9 (1.0 til) X 10"^ (16) 

with the non-zero best fit value. 

The reconstructed matrices for the normal, inverted mass hierarchy and the degenerate 
mass spectrum are given in Table 1, where we show the absolute values of the matrix 
elements. For each case we present (i) the matrix for the best fit values of the parameters, 
(ii) the intervals which correspond to the la experimental uncertainties for zero CP-phases 
and the intervals when also the CP- violating phases vary in whole possible range: Xi — O — ir, 

The following comments are in order. 

1) In the case of normal mass hierarchy we take mi = 0, so that m2 = y^Am^i and 
777,3 — y'Amli. Notice that the 1-3 element of the matrix is much smaller than the 1-2 
element and elements of the dominant block are different. This is the consequence of non- 
zero (though statistically insignificant) 1-3 mixing and shift of the 2-3 mixing from maximal 
one. The la experimental uncertainties lead to SrUee ~ 2.5 meV, 5me^ ~ 5mer ~ Sm^r ~ 5 
meV, (5m^^ ~ drrirr ~ 10 meV. The effect of CP-phases is sub-leading. Variations of phases 
A3, S result in similar size of the intervals and both effects double the indicated uncertainties. 

Notice that there are strong correlations between elements, so that their values can not 
be taken from the indicated intervals independently. 

2) For the inverted hierarchy we take 7773 = 0, mi f» m2 ~ -y/jAm|7|- We choose 
Ai = 5 = (A3 is irrelevant) for the best fit data analysis. The la experimental uncertainties 
lead to Sm ~ (8 — 10) meV for all matrix elements but 5m^,- ~ 4 meV. The effect of 
phase variations is much stronger. Now the phases affect the dominant block and therefore 
^777 ~ 777. 
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3) Degenerate spectrum. For illustration we take mi ^ m2 ~ m^ = m^ = 0.2 eV. The 
experimental errors produce very small effect: Sm ~ 1 meV. In contrast, effects of phase 
variations are very strong, 5m ~ mo- Here also variations (values) of different elements are 
strongly correlated. 

2.3 Induced matrix and the dominant structures 

The mass and mixing patterns in quark and lepton sectors are strongly different. The 
difference can (at least partially) originate from the active-sterile mixing which is absent in 
the quark sector. The shortcoming of this proposal is the coincidence problem: two different 
contributions to the mass matrix, active, ma, and induced, mi, are of the same order or 
within 1-2 orders of magnitude in spite of the fact that they have different, and at the first 
sight, unrelated origins. The only argument in favor is that this will be not the only case - 
we meet the coincidence problem in other areas too. Another possibility is that mi ^ ma. 

Let us consider first the case of one sterile neutrino. Due to the singular (rank 1) 
character of the induced matrix it can not reproduce the dominant structures of the neutrino 
mass matrix in the case of the degenerate mass spectrum and inverted mass hierarchy. In 
the former case Detm ml, where mo is the scale of neutrino mass. In the latter - there 
are two dominant eigenvalues and the determinant of the 1-2 submatrix is non-zero. 

Essentially this means that the matrix induced by one sterile neutrino can be the origin 
of the dominant block only in the case of normal mass hierarchy. This can explain the large 
(maximal) 2-3 mixing. Suppose 



m^s = rrirs = mo, 



mes -C mo. 



(17) 



then 




















1 



1 



(18) 



ms 



\ 







1 




From (18) we obtain 2ml/ ms — V^^ai 





(19) 
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Parameters of the 1-2 block (mixing and mass) should be given by the original active 
mass matrix. For the best fit values of the oscillation parameters (here we take sin^ ^13 = 0) 
and the induced matrix (18) the active neutrino mass matrix should be of the form 

\ 

(20) 



U^, U!t m2 Ult^ 



'23 



m2 



V2 

^2 



V 



sc 
x/2 

ci 
2 
.2 



where m2 = diag(0, m2, 0), s = sin ^12, c = cos^i2, 11^2^ is the 1-2 rotation matrix on the 
solar mixing angle, f/23 the matrix of maximal 2-3 mixing. We assumed that mi ~ 0. 
Notice that all elements of this matrix are nearly equal being in the range 0.31 - 0.35. 

A non-zero mi would be equivalent of adding to (20) the matrix mil proportional to 
the unit matrix and substituting m2 (m2 — mi), if there is no CP- violating phases. That 
does not change our conclusion provided that mi ^ m2. 

The induced matrix can be chosen in such a way that the active one has hierarchical 
structure with small mixings similar to the charged fermion mass matrices. (Partly this case 
has been studied in [6]). Consider non- universal active-sterile coupling of the type 



rrias = ("^3 - mi)ms {a, a, b), a <^ a b. 



Then the required mass matrix of active neutrinos can be written as 



(21) 



nia = mil -|- (m3 — mi) 



V 



es^ — o? 



e-^-aa 



-e^ - ab 



1 



ab 



1-b^ + e 



(22) 



2 / 



where e = m2/ (m3 — mi). For a — 0.95 and b — 1.50 we obtain hierarchical 2-3 block with 
("^0)22 '■ ("^a)23 '■ ("^0)33 ~ 0.1 : 0.4 : 1.2. Other elements are of the order 0.02 - 0.04 if 
a ~ 0.01. Certain hierarchy among those elements can be obtained for non-zero 1-3 mixing 
and some deviation of the 2-3 mixing from maximal. The overall spread of the values of 
elements by 2 orders of magnitude can be easily obtained. 

Another scenario (essentially considered in [6]) is that 6 ^ a, so that the induced matrix 
corrects the 33 element only. Taking, e.g., rriss — sin^^^^m^ ~ 300 meV and the final 



matrix as in the Table 1, we obtain that the original active mass matrix should be like 
that for normal hierarchy (Table 1) but with 777.33 ^ 300 meV, that is, strongly hierarchical. 



2.4 Induced matrix and mass hierarchy 

An interesting possibility is that the induced matrix can switch the mass hierarchy from 
normal to inverted and vice versa. Indeed, 



V2 



m. 



|A777i2 



D, 



(23) 



where the induced term, D, is close to the "democratic" matrix with all elements being 
nearly 1. 

2.5 Induced matrix and tri-bimaximal mixing 

The mass matrix which generates the tri-bimaximal mixing in the normal mass hierarchy 

case can be presented as 



111 = 



"^32! 







V 



+ 



A777^ 



^ 1 1^ 



21 



111 

1 1 1 



(24) 





1 -1 
-1 1 

It is the sum of two singular matrices. The sub-dominant (second) matrix can be induced 
by universal mixing with the sterile component 



maS = "7o(l, 1, 1). 



Then according to (24) 



7775 sm Us 



li 



3 meV. 



(25) 



(26) 



In fact, both matrices in (24) can be induced by the active-sterile mixing, if the second 
sterile neutrino is introduced with mixing elements 

777:^5 = 7770(0,1,-1). (27) 

In this case the original active neutrino masses should be very small, e.g., of the order of 
the Planck mass suppressed scale: ~ v^^/Mpi. 
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2.6 Induced matrix and QLC 

If not accidental, the quark-lepton complementarity (QLC) relation [14] can imply that 
(i) there is some structure in the lepton sector which generates the bi-maximal mixing 
f/bm = f^23 ^125 (ii) there is the quark-lepton symmetry which "propagates" the CKM- 
type rotations to the lepton sector. 

Let us consider a possibility that the induced matrix is responsible for the bi-maximal 
mixing: 

mi = nib™ = U^^m^'^^U^^: (28) 

whereas the charged lepton mass matrix produces the CKM-type rotation. The original 
active neutrino matrix should then give very small contribution. 

Clearly this scenario can not be realized with only one sterile neutrino: Taking m*^'^^ = 
diag(0, 0, ms) we find that mj = U^iUsU^^ with a single maximal mixing. With two 
sterile neutrinos the matrix nibm which generates the bi-maximal mixing can be reproduced 
precisely. As an example one can take nias — mo{a,b, —b) and m'^g — mo(0, a;,y), where 
X — \/d-\- Q.ba? — 6^, y — {d — 0.5a^ -|- 6^)/x, and a, 6, d are free parameters. 

2.7 Small and negligible induced contribution 

The induced matrix becomes irrelevant if 

« lm.hr (29) 

ms 

Let us find conditions under which the effects of sterile neutrino are below the present 
1(7 spread of the matrix elements. Uncertainties depend on the type of mass spectra and 
are different for different elements. 

For the normal mass hierarchy the effects are below la if 

sin^ 9es < 2 meV, sin^ ms, sin^ 9t-s < 5 meV, (30) 
For II and r we have taken the uncertainty of m^,- which is the smallest one. 
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For the inverted hierarchy the induced contributions are below la uncertainties if 

sin^ 9es ms <S meV, sin^ 6*^5 777.5 < 4 meV. (31) 

Larger effects of sterile neutrinos, sin^ ^57775 < 20 meV, can be mimicked by the phase 
variations. 

In the case of degenerate spectrum la experimental uncertainties restrict 

sin^ Oes ms <1 meV. (32) 

The phase change produces the same effect as sterile neutrinos with sin^ 9as ~ 200 mcV. 

Sterile neutrinos can be responsible for fine structures of the mass matrix. The smallest 
(observable) structure is related to the solar mass split. In the case of inverted hierarchy 
this would correspond to the contribution sin^^g^mg ~ 1.5 meV, and for the degenerate 
mass spectrum we obtain the smallest quantity: sin^ OesiTis ~ 0.4 meV. 

So, we can identify three benchmarks: 

1) For sterile neutrinos with mixings and masses smaller than 

sin^ 6as ms — 1 meV (33) 

the effects are below the present la experimental uncertainties for the hierarchical spectra. 
Still these neutrinos can influence the sub-leading structures in the case of the degenerate 
spectrum. 

2) Sterile neutrinos with 

sin^ e^s m5 = 3 meV (34) 

can generate the sub-leading structures in the case of normal mass hierarchy. 

3) Sterile neutrinos with 

sin^ 0^s ms = (20 - 30) meV (35) 

can generate dominant structures in the case of normal and inverted hierarchies. For larger 
masses and mixings they can lead to dominant structures of the degenerate spectrum. 

In Figs. 1, 2 and 3 we show the lines of constant induced mass sin^ O^s ms — const., in 
the plane sin^ 9as and 7775 which correspond to the values in (33, 34, 35). 
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3 Bounds on the active-sterile mixing 



In this section we describe the direct mixing effects of S and bounds on its masses and 
mixings. 

3.1 Production, thermalization, decay 

The most stringent bounds follow from astrophysics and cosmology. Sterile neutrinos can be 
produced in the Early Universe non-thermally through their mixing with active neutrinos, 
affecting primordial nucleosynthesis, cosmic microwave background radiation as well as the 
growth of cosmological structures. 

We assume that no primordial density of sterile neutrinos existed, and all sterile neutrinos 
where produced in the Early Universe due to mixing with active neutrinos and oscillations [9, 
10]. 

The "thermahzation" lines and the "decay" lines in the ms — sin^ 9s plane (Figs. 1-3) 
allow to understand various bounds. 

1) The thermalization lines give the lower bounds of the ms — sin^ 6s region where 
the sterile neutrinos are thermalizcd before the primordial nucleosynthesis. According to 
Ref. [15, 16] the lines are slightly different for mixing with electron and non-electron neu- 
trinos: 

ms sin^ 9s — 0.6 meV, for Ue, (36) 

ms sin^ 9s — 0.4 meV, for v^^Vr- (37) 

Notice that the thermalization lines have the same functional dependence as the isolines of 
induced mass. Furthermore, the lines are below the benchmarks obtained in (33, 34, 35). 
That is, sterile neutrinos with parameters which give significant induced contribution were 
thermalized in the Early Universe. 

2) We confront the life time of neutrinos with the time of recombination. Tree ~ 10^^ 
s, and the age of the Universe, ru = 4 x 10-*^^ s. The decay rate of S strongly depends 
on ms- For small ms the main channel is S ^ ?>v and then the following channels open: 
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5— + / + 5— s>z/ + g + g, etc.. The decay rate can be written as 



1 



's, 



(38) 



where and are the muon mass and decay rate, correspondingly, and K{ms) is the num- 
ber of decay channels for a given value of uis. So, for the non-relativistic sterile neutrinos, 
the lines of constant decay rate are given by 



In Figs. 1, 2 and 3 we show two decay isolines corresponding to ts = Tree and Tg = tjj. 
3.2 LSS formation bound 

If Ts > Tjj, the sterile neutrinos contribute to the dark matter in the Universe. Analysis of 
the Large Scale Structure (LSS) of the Universe gives the bound on the total energy density 
in sterile neutrinos, ps, as function of its mass: 



Here Pcr is the critical energy density and h is the Hubble constant"^. We use for ujsims) 
the results of the analysis in [17] and [18]. 

For ms < 100 eV, S compose the hot dark matter component and the bound on us is 
stronger. According to [17] for ms < 30 eV one has (Vs < 0.005 at 95% CL. The bound 
weakens with the increase of ms, as neutrinos become colder: a;^ < 0.02 at ms = 100 eV, 
and Us < 0.12 at ms = 300 eV. In the interval ms = (0.25 - 30) eV the 95% CL bound 
from [18] can be parametrized as 



n{ms) m% sin^ 6s = const. 



(39) 



ujs = — < uJs{ms). 

Per 



(40) 




(41) 



^In this paper in numerical estimations we use h=0.7. 
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For ms = (30 - 300) eV ^: 



ujs < 0.01 (in (^) - 2.55 



2.2 



V Vev; 



- 1.0 



(42) 



and ojs < 0.12 for ms > 300 eV. 

Calculations of the energy density a;^ have been updated recently [19] with inclusion of a 
number of additional effects, in particular, effects of the quark-hadron transition, modifica- 
tion of the finite temperature effective mass of the active neutrinos, etc.. For the temperature 
of the quarks-hadron transition Tqcd — 170 MeV, the relation between cus, ms and 9s can 
be parametrized as [19] 



Notice that for ms < 3 keV it deviates from the dependence previously found in [20]. 

Plugging the hmits (41) and (42) in relation (44) we find the upper bound on sin^ 9s as 
a function of ms (see region LSS in Figs. 1-3). The bound is absent for ms < 0.25 eV [21]. 
It docs not depend on the flavor of active neutrino to which sterile neutrino mixes. 

The bound is valid for the region of parameters below the isoline Ts = Tu which corre- 
sponds to ms < 0.5 MeV and sin^ 9s > 7 ■ 10~^^. However for ms > 5 keV the stronger 
bound follows from the diffuse background radiation. 

It has been shown in a recent analysis [22] that sterile neutrinos with mass ms < 14 
kcV arc excluded at 95% CL as the dark matter particles responsible for the LSS formation. 
This implies somehow stronger bound on uJs than the one we use in our estimations and 
consequently stronger bound on the S parameters. It is argued in that the bound can be 
relaxed if new particle decays in the epoch between the decouphng of S and BBN increase the 
*A more accurate parametrization in this range is: 



0Js<Q- 10-3 + 4.1 • lO"'^ (ms/eV) + 2.5 • 10"'^ {ms/eVf + 1.7 • 10"^ {ms/eVf - 4.5 • 10"" {nis/eVf 




(43) 



or 




(44) 
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entropy, thus diluting concentration of S and reducing their relative temperature [23] . This 
may be the case for the model proposed in [24, 25] where two heavy RH have masses above 
1 GeV. (Explanation of LSS by S would probably require that S are produced at some high 
energy scales by processes that are not related to active-sterile neutrino oscillations [23] .) 

3.3 Limit from cosmic X-ray radiation 

Due to mixing via the loop diagrams sterile neutrinos decay into an active neutrino and a 
photon: S with Ri rns/2 and at the rate ~ a^su- Therefore one expects to 

detect the photon emission line when looking at big concentrations of dark matter such as 
galaxy clusters. 

The analyzes of the X-ray emission from the Virgo cluster [26, 27] give limits on the 
decay and therefore on the mixing of sterile neutrino: sin^^g < 2.6 x 10~^ (m^/keV)"^. 
(This parametrization is valid for ms = (1 — 10) keV.) It is argued in [28] that the bound 
from Virgo is about 1-2 orders of magnitude weaker especially for low masses, ms < (5 — 8) 
MeV. However for ms — 10 keV, relevant for this analysis (where X-ray bounds start to 
dominate over the LSS one), the bounds are comparable: sin^ ^5 < 2.6 x 10~^° [26] and 
sin^ 9s < 5.0 x 10"^° [28]. At the same time stronger hmit has been obtained from analysis 
of Coma cluster [28]: 



The 7-flux from all possible sources could accumulate over the history of the Universe 
and be seen as a Diffuse Extragalactic Background RAdiation (DEBRA). Apparently the 
flux from the radiative S decay should be smaller than the observed flux. In the range 
ms = (1 — 100) keV the DEBRA limit can be parametrized as [29] 



This bound based on the data collected from the whole sky is weaker than the one given 
in (45). However, it does not depend on assumptions concerning clustering and therefore 
is considered to be more robust [29]. These limits are valid below the recombination hue 

TS = Tree- 




(45) 




(46) 
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The X-ray exclusion region shown in Figs. 1-3 has been obtained from the analysis of 
Coma cluster data as well as data on diffuse X-ray background (DEBRA) from HEAO-1 
and XMM-Newton missions in [29]. We have also included the very recent bound from the 
diffuse X-ray spectrum of the Andromeda galaxy [30] . This is the most stringent limit in 
the range ms = (1 — 24) keV. For small masses, ms < 1 keV, we use the hmit from [31]. 

Note in Figs. 1-3 that the X-ray data reduces the parameter space allowed for Warm 
Dark Matter to a very small region: rus = (1.7 — 3.5) keV for sin^ 9s ~ 10^^ — 10~^. Also 
most of the parameter region that can explain the origin of pulsar velocities is ruled out, 
remaining only a corner around ms ~ 2 — 4 keV. 

It was proposed in Ref. [26] to observe clusters of galaxies with Chandra and XMM- 
Newton observatories, in their high sensitivity range for X-ray photon detection of (1 - 
10) keV. That will allow one to set the limit in the range 10"^"^ < sin^2^5 < 10~^ for 
ms = (0.6 - 40) keV. 



If the sterile neutrinos decay producing light neutrinos between the active neutrino decou- 
pling time, Tu ^ 1 s, and the photon decoupling time. Tree, this would increase the energy 
density of relativistic particles at t < Tree. The density is described by the effective num- 
ber of neutrinos, N^y. That, in turn, affects the CMB angular power spectrum (acoustic 
peaks). The bounds on N^, from observations were substantially improved during the last 5 
years: < 13 at 95% CL (BOOMERanG/MAXIMA) [32], < 8.3 at 90% CL (WMAP 
data) [33], < 6.8 - 7.1 at 95% CL (WMAP and LSS) [34, 35, 36], < 5.4 at 95% CL 
(WMAP, LSS and type la supernova data) [37] . 

With the help of considerations in Ref. [38] and [20], we get the "CMB" limit on the 
sterile neutrino parameters as a function of A^i^, which can be parametrized as 



Combined analysis of the cosmological data on LSS, supernovas and the CMB including 3 



3.4 CMB-bound 




(47) 
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years result from WMAP [39] allows to put the bound Ny < 3.74 at 95% CL. Using this 
results and Eq. (47) we find the bound shown in Figs. 1-3. This limit is not valid above the 
fine Ts < T„ ^ 1 s, which for sin^ 9s — 10~^^ corresponds to ms — (400 — 500) MeV. 

In the future, the PLANCK mission [40, 41] will allow to strengthen the bound down to 
< 3.2, while according to Ref. [41], the CMBPOL mission can achieve < 3.05. This 
will further expand the excluded region, in particular to larger values of ms- 

3.5 BBN bound 

Apparently for one sterile neutrino the limit on the effective number of additional degrees of 
freedom during the epoch of Big Bang Nucleosynthesis (BBN), AA^,^ = 1, does not provide 
any bound since at most the equilibrium concentration of S can be produced in the scenarios 
under consideration. The limit AA^,^ < 1 becomes relevant in the case of more than two 
sterile states. 

For the low values of ms we use the limits from [42] which for AA^i^ = 1 can be 
parametrized as 

7715 sin^ 6*5 = 1.4 meV, for u^, (48) 

sin^ 6's' = 1.0 meV, for f^iVr- (49) 

For the high masses, ms = 10 — 200 MeV, the exclusion region has the shape of a 
parallelogram [43] . The right boundary is basically given by 

which is the condition for the heavy neutrinos to be relativistic at decoupling so their 
number density is not Boltzmann suppressed and they can have an impact on BBN. The 
left boundary, 

can be simply understood as the condition that the energy density of the sterile neutrino 
is smaller (including also an entropy dilution factor 5) than the energy density of one light 
neutrino species at BBN. 
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3.6 Supernova neutrino bound 

Two different bounds follow from observation of the antineutrino signal from SN1987A [44]. 
One is from Pg-disappearance and the other from star cooling. 

1) The resonance conversion i?e ^ S can occur in the central regions of the star due to 
change of sign of the matter potential. That happens for the active-sterile system when the 
relative electron number density is Ye ~ 1/3. For the mass range ms ~ (1 — 100) eV the 
adiabaticity condition can be fulfilled if sin^^g > 10~^. The adiabatic conversion leads to 
strong suppression of the P^,- flux [3]. So, the observation of the Pg-signal from SN1987A 
gives the bound on the oscillation parameters. In Figs. 1-2 we show the updated results 
obtained in Ref. [4]. 

2) For large masses, ms, the {ua — > S) oscillations as well as scattering lead to production 
of sterile neutrinos in the core of the collapsing star. If these neutrinos escape the core, they 
can lead to substantial energy loss. This, in turn, will shorten the neutrino burst and the 
energy released in z/g will be smaller [31, 45]. 

Normal duration of the SN1987A neutrino burst excludes significant cooling effect which 
puts the bound on S parameters [43] shown in Figs. 1-3. Notice that the upper bound, 
sin^ ^5 < lO"^*^, is shghtly different for and and u^-. If mixing is sufficiently large, 
the sterile neutrinos will be trapped inside the core, and the energy-loss argument is not 
applicable. This gives the lower bound of the excluded region sin^ ^5 < 10~^. The mass 
range is restricted by the condition that S is produced inside the core. 

Future detection of high statistics neutrino signal from a Galactic supernova will allow 
to put stronger bounds. Indeed, another resonance in Ue — S can occur in the outer regions 
of the star with normal chemical composition. The adiabaticity condition can be written as 



which in fact is stronger than the cosmological bound in the range (0.01 - 0.1) keV. The 
adiabatic — S conversion leads to the disappearance of the i/g-neutronization peak, and 
the modification of the signal during the cooling stage. In particular, in the case of normal 
mass hierarchy the electron neutrino fiux at the Earth will be -F(t'e) = F^ii^n) in the case 




(50) 



18 



of large 1-3 mixing and -F(fe) = cos^ ^i2-F°(t'^) for very small 1-3 mixing [3]. 

In the eV mass range, sterile neutrinos can drive non negligible i^e/^n ~^ ^ei^r 

conversions at TeV energies by the MSW effect that can be constrained by future IceCube 
data [46]. 



3.7 Laboratory bounds 

The laboratory bounds are typically much weaker than the astrophysical and cosmological 
ones. They are, however, more robust and turn out to be the main bounds if, for some 
reason (see sec. 4.2), the cosmological and astrophysical limits become inapplicable. 

1). The Ou^P- decay. Introduction of active-sterile mixing (described by nias) does 
not change the rriee element of the whole 4x4 mass matrix (4). Therefore for light S: 
ms <S i/fN, where tat is the typical size of nuclei, the rate of Oi//3/3-decay in the lowest 
approximation is determined by iriee and the effect of S is strongly suppressed. 

However, the effect of S on the 0z//3/3-dccay increases with ms- Let us consider this in 
more details. The QuPf] decay amplitude, A^^\ has two contributions associated to S: i) 
from the induced mass, {mi)ee, and exchange of light neutrinos, and ii) from exchange of ^S" 
and its mixing with i/e- So, the amplitude of the decay can be written as 

— mf, — rrig 

where g ~ l/r^r is the effective momentum of the exchanged neutrino. As it can be inferred 
from the calculation presented in [47], q pa 100 MeV. In the denominator of (51) the pa- 
rameter nil, is the effective mass of light neutrinos and since rrii, <S ms it can be neglected. 
Taking into account that the induced mass {mi)ee — —ms sin^ O^s we can write, according 
to Eq. (51)), the total contribution of S to the effective Majorana mass as 



m^f ) = ms sin^ O^s 



1- 



q^ — ms 



(52) 



If ms is negligible the two contributions cancel each other leading to zero effect of S. If 
m| <^ g^, we find from (52) 

m(f)=sin^e,5^. (53) 
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So, for light S, ms <^ '^/^n, ~^ ^'^d the total effective Majorana mass is determined by 
the ee-elements of the active neutrino mass matrix (1113)^^, in agreement with consideration 
in terms of 4 x 4 matrix. 

For m| ^ g^, the second term in (52) can be neglected and the Majorana mass is given 
by the induced contribution (as in the usual see-saw mechanism): 

mfj = sin^ e^sms ■ (54) 

In Fig. 1 we show the excluded region of the S parameters which corresponds to the upper 
bound mee < 0.5 eV obtained from studies of Oi//3/3-decay of ^^Ge [48]. We assume that 
there is no cancellation between S contribution and (ma)gg. In general, cancellation reduces 
the excluded region. However, in our context, the cancellation can not be strong. Indeed, 
according to our assumption (Sec. 2.1) the original active neutrino mass terms are below 1 
eV. Therefore maximal contribution to the effective Majorana mass of electron neutrino is 
about 1 eV, and consequently, only the order of 1 eV contribution from new neutrino state 
can be cancelled. That corresponds to the long-dashed line in Fig 1. 

As follows from Fig. 1 the double beta decay limit becomes relevant for high masses 
ms > 100 MeV. 

Future neutrinoless double beta decay experiments [49] with sensitivity down to < 
(0.01 — 0.03) eV will improve the bound on sin^ 9s by a factor 10-30. 

2) The 13-decay. The region ms ~ (0.1 - 10^) keV and sin^ Og < lO'^ - lO'^ for i^s - Ve 
mixing (Fig. 1) is excluded by the negative results of searches of kinks in the energy spectra 
of nuclear /3-decays [50] . 

3) Meson decays can have contributions from sterile neutrinos which can modify the 
energy spectra of their decay products. In the us — i^e channel, the best limit comes from 
precision measurements of the energy spectrum of e"*" in the decay tt"*" e^Vf. [51]: sin^ 6s < 
10~^ for ms — (50 — 130) MeV. In the us — v^j, channel, the bounds come from studies of the 
spectra of 77=*= and decays at accelerators and in the atmosphere. The 90% CL excluded 
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region (labeled as "Decays" in Fig. 2) is taken from Ref. [52]. It is the only bound which 
can compete with the cosmological and astrophysical limits in the range of large masses 
ms ~ (0.03 - 0.3) GeV. 

In the near future, the studies of pion, muon and kaon decays at MiniBooNE [53], 
MINOS [54] and K2K [55] experiments are expected to reach sensitivities of a fewxlO"^ or 
less for the mixing if nis ~ 100 MeV [52]. 

Oscillations of active to sterile neutrinos lead to the suppression of the neutral current 
interaction rate. In this connection, it has been proposed to study interactions of neutrinos 
from pion decays at rest with the superallowed neutral current reaction 1/^ ^^C — > ^^C* [56]. 
The sensitivity of such an experiment can reach sin^ 9s ~ 10~^ for ms ~ 1 eV. 

4) The atmospheric neutrinos and K2K. The combined analysis of the SuperKamiokande [57] 
and MACRO [58] results as well as the data from accelerator experiment K2K leads to the 
following 90% CL. bounds [4]: for the — vs channel sm^Os < 10~^, in the interval 
777,5 = (0.01 — 0.1) eV (see Fig. 1); for the i/^ — us channel sin^ ^5 < 0.04 in the range 
'^s = (0.1 — 1) eV (see Fig. 2); for the u^- — z/5 mixing sin^ < 0.1. 

If the cosmological bound is invalid, the atmospheric bound becomes also relevant for 
larger masses, 7775 = (1 — 10) eV. 

Sterile neutrino oscillations for 7775 ~ few eV have the oscillation length comparable to 
the Earth's radius at ~ 1 TeV. In the future, IceCube [59] detection of atmospheric 
neutrinos can probably bring the limit on i/^ — i>s mixing down to ~ 10~^ for ms ~ few 
eV [46]. 

It has been suggested in Ref. [60] that showers generated by ultra-high energy sterile 
neutrinos {Ej, ~ 10^ — 10^^ GeV) with sin^ 9s ~ 0.01 — 0.1, may be distinguished from the 
ones generated by active neutrinos in experiments such as EUSO [61] or OWL [62], using 
the Earth as a filter, as proposed in Ref. [63] for neutralino showers. The mass range of ms 
that can be tested by these experiments is determined by the mechanism of their production. 

5) Reactor neutrino experiments can set limits on the mixing Pg — P5 by comparing 
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expected with observed Ve flux. Bugey [64], CHOOZ [65] and Palo Verde [66] provide the 
main bound on — us mixing, sin^6's' < 0.01 (90% CL.), in the range ms = (0.3 — 1) eV. 

The forthcoming reactor experiment Double-CHOOZ and the proposed projects Kaska, 
Braidwood and Angra [67], may be able to reach sin^ ^5 ~ 5 • 10~^ for ms — (0.03 — 0.7) 
eV [68]. 

6) Accelerator bounds. The accelerator oscillation experiments CDHS [69], CCFR [70], 
NOMAD [71] and CHORUS [72] as well as LSND [8], KARMEN [73], provide with the 
bound on 1/^ — S mixing in the range low mass range, ms > (1 — 100) eV. In particular, 
sin^^s < 7 • 10-3 (90% CL ) for ms > 10 eV (see region labeled "Beam" in Fig.2). 

In Figs. 1 we show the forbidden region of the parameters (labeled "Reac.+Beam") from 
the combined analysis of these accelerator beam and reactor experiments Ref. [4]. This 
result excludes the induced mass as the origin of the dominant structure in the low ms 
domain. 

In the high mass range, ms > 0.1 GeV mixing with a heavy S leads to an effective 
violation of lepton universality and appearance of the flavor changing neutral currents. The 
non-observation of these effects in experimental data permitted the authors of Ref. [74] to 
set the 90% CL hmits, shown in Figs. 1-3: sin^ 9s < 7.1 • 10"^ {us — Ve), vahd in the range 
ms > 0.14 GeV (m^r), sin^ 6*5 < 1.4 • lO^^ {vs - ^ix), valid in the range ms > 1.115 GeV 
(ttia), and sin^^5 < 1.7 • 10"^ (1^5 — i/^), valid in the range ms > 1.777 GeV (m^). Also 
searches for a singlet neutral heavy lepton at LEP [75], have allowed to exclude sin^ 9s down 
to 10"^ — 10"^ depending on the channel, in the range 0.4 < ms/GeV < 90 at 95% CL. 

Accelerator experiments sensitive to oscillations with Am^ > Am|2 will be able to set 
stringent bounds in the low mass range. The MINOS detector [54] can use the neutral 
current to charged current ratio to probe — us and Ur — i^s mixing down to sin^ 9s ~ 
10"^ [4]. Sensitivity to v^ — i's mixing at the level sin^ ^5 ~ 10 ^ for ms = (0.03 — 0.3) eV can 
be reached by T2K [76], NOi/A [77] and future z/-factories looking for i/^, disappearance. 

Forthcoming pp and planned e^e~ collider experiments can expand the the excluded 
region of S parameters to larger masses. According to Ref. [78], the search for same-sign 
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dilepton production mediated by a sterile neutrino S in pp ^ l+l +X with /, /' = e, /i, r, may 
be used to constrain sin^ ^5 < 10~^ for ms = (0.1 — 2) TeV at the Large Hadron Colhder 
at CERN, provided a large integrated luminosity becomes available. It has been shown [79] 
that at an International Linear Collider with a center of mass energy ^/s — 500 GeV, one 
can look for single heavy S production through e+e~ Sv ^ IWv with I — e,jji, r. This 
has a sensitivity down to sin^ 9s < 7 ■ 10"^ for ms = (200 — 400) GeV. The same reaction 
at a future Compact Linear Collider operating with ^/s = 3 TeV could be used to limit 
sin^ ^5 < (2 - 6) • 10-3 for = 1 - 2 TeV [80]. 

Notice however, that these high mass bounds are essentially irrelevant in our context. 
Indeed, for ms > 1 GeV and sin^^^ > lO^^, the induced mass is m/ > 1 MeV. This 
means that the elements of the original active neutrino mass matrix should also be large 
ma > 1 MeV, and it should be extremely strong cancellation of the original and induced 
contributions to obtain phenomenologically acceptable masses of light neutrinos: ((ma)y — 
(mi)y)/((ma)y < 10-^. 

Inversely, discovery of the S with relatively large mixing in high energy collisions will 
testify against the approach developed in this paper. 

4 Induced mass versus direct mixing effects 
4.1 Bounds on induced mass 

Confronting the lines of constant induced masses (33), (34) and (35) with cosmological 
astrophysical and laboratory bounds in the sin^ 9s — ms plane we can conclude on the 
relative importance of the direct and indirect effects. 

There are two regions of parameters in the sin^ 9s — ms plane, where the induced masses 
are more important than the effects of direct mixing. That is, in these regions substantial 
induced masses for the active neutrino mass matrix are not excluded by the existing bounds. 

1) High mass region: ms > 300 MeV and sin^^s < 10~^. This region is restricted 
essentially by the CMB bound, meson decays and SN1987A cooling. Future measurements 
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can probably improve the bounds by about 1 order of magnitude from below. 

Only for the Vf. — vs channel, the neutrinoless double beta decay can probe the whole 
high mass region. Present bounds correspond essentially to the dominant contribution in 
the case of the degenerate spectrum. 

For the other mixing channels, it is the induced mass which gives the bound on the 
parameters of the sterile neutrinos. Indeed, assuming that there is no strong cancellation 
of elements of the original active neutrino mass matrix nia and mi, and taking the largest 
elements of matrices in Table 1, we can write the bound 

sin^ Qs-ms < m^^P - (0.5 - 1) eV, (55) 

which is clearly comparable to the neutrinoless double beta decay bound but now valid for 
all channels. This can be viewed as the forth benchmark line shown in Fig. 1-3. 

Sterile neutrinos in this range can play some role in leptogenesis and the generation of 
the baryon asymmetry in the Universe [81, 24]. 

Clearly contribution from this region is out of our control and this creates ambiguity in 
the implications of the mass and mixing results. 

It is beyond the scope of this paper to discuss the possible origins of S with such a 
small mixing. Still one possibihty looks rather interesting: if ms — (10^ — 10^) GeV (region 
where one may expect singlets related, e.g., to SUSY breaking), the required mixing is 
sin^ 9s = 10^^^ — 10^^^. The latter can be related to the existence of a new "intermediate" 
scale M = ms/sinOs ~ (lO^o - lO^^) GeV. 

If the interpretation of the LSND result in terms of oscillations in (3 + 1) scheme is 
confirmed, that would imply existence of the sterile neutrino (s) with mass (0.5 — 5) eV with 
mixing parameters sin^ 9s ~ 0.02 (see [82] for recent analysis). The corresponding induced 
mass equals m/ = (10 — 100) mcV and therefore the effects of LSND neutrino on the active 
neutrino mass matrix is strong and can not be considered as small perturbation. 

2) Low mass window: nis ~ (0.1 — 0.3) eV and sin^^5 = 10~^ — 10~^. This window 
is essentially closed for all the channels if one takes the BBN bound AA^,^ < 1. In this 
case, one has the bound on the induced mass m/ < 1 meV. That can produce some effect 
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on observablcs in the case of degenerate or inverted hierarchy spectrum. From above this 
region is restricted by the LSS bound on neutrino mass. 

If AA^,^ = 1 is allowed, there is no BBN bound. Bounds from other effects strongly 
depend on flavor. For the — vs channel, the reactor and atmospheric neutrino bounds es- 
sentially exclude the dominant contribution from m/ but still allow for sub-dominant effects. 
For the z/^ — vs channel, the bound is given essentially by the atmospheric neutrinos and 
larger region of mixings is allowed. In particular, for ms ~ 0.25 eV the mixing sin^ Os = 0.04 
is not excluded leading to m/ ~ 10 meV. The latter is close to the dominant contribution. 
For the 1/^ — channel, the atmospheric bound is weaker and dominant contributions from 
mi are allowed: mi ~ (30 — 250) meV. 

In the rest of the ms region, ms = (10~^ — 10^) keV, effects of the direct mixing dominate 
over the induced matrix effects. In this range the induced masses 



can produce only very small corrections to the active neutrino mass matrix. In the interval 
ms = (1 — 10^) keV the bound is even stronger: 



If however the cosmological/astrophysical bounds are absent for some reason, a large 
range of parameters becomes allowed and the induced masses can reproduce the dominant 
structures of the active neutrino mass matrix in the whole range of ms- 

In Fig. 1 - 3 we show also regions of parameters which correspond to certain positive 
indications of the existence of sterile neutrinos: (1) overlapping regions of the warm dark 
matter and pulsar kick, (2) the LSND spot. 

In the scenario [24] two neutrinos have masses ms ^ 1 GeV and their mixing is respon- 
sible for the mass matrix of hght active neutrinos: nia = m/. The third sterile neutrino 
with mass ms ~ 1 keV and mixing sin^ 6s ~ 10~^ — 10~^ can contribute substantially to 
warm dark matter (WDM) and explain pulsar kicks. 



mi<A- 10"^ meV, 



(56) 



mi < 10 ^ meV. 



(57) 
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Apparently the WDM and pulsar kick regions are far below the benchmark lines and 
therefore the corresponding induced masses are negligible. Furthermore, the regions are 
disfavored by the cosmological and astrophysical observations. The WDM scenario can 
be recovered if the mixing is smaller then that indicated in the plot and some additional 
mechanism of generation of sterile neutrinos exists apart from mixing with active neutrinos. 

In contrast, the LSND spot is in the range where the induced masses are of the order 
of dominant mass structures. So mixing with sterile neutrino can not be considered as 
small perturbation of the original active neutrino structure. The LSND spot is essentially 
excluded by the cosmological data unless some new physics is added. 

4.2 Avoiding bounds 

Let us consider various possibilities which allow one to circumvent the bounds obtained in 
the previous section and therefore to open a possibility for strong effects of the induced 
matrix even for low m^. 

1) It has been shown [60] that in a cosmological scenario with low reheating temperature 
at the end of inflation, Tr << 100 MeV, experimental bounds on the active-sterile neutrino 
mixing are relaxed. In particular, the sterile neutrino required to explain the LSND result 
is allowed in this scenario. Also the influence of a non-negligible primordial lepton number 
asymmetry has impact on the sterile production rate since in this case neutrinos are produced 
resonantly with a non-thermal spectrum [83, 20, 84] . According to Ref . [20] , for L = 10~^ 
the LSS hmit presented in Figs. 1-3 will change very little, while for L — 0.1 the hmit on 
sin^ 9s can increase by 1 or 2 orders of magnitude depending on ms- 

2) Another possibility can be related to the origin of the mass of S itself. Recall that 
S may not be related to right handed neutrinos and usual family structure. Let us assume 
that S has a "soft mass" generated by the medium dependent VEV of some new scalar 
field A: ms — ^{A). The VEV can be proportional to the number density of the active 
neutrinos, Ui,: {A) oc n^, [85]. In this case 

ms^m%{l + z)\ (58) 
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where ni^g is the mass in the present epoch. If the mixing mass is generated by the usual 
Higgs VEV, which does not change with time, we find that the induced contribution to the 
fight mass and the active-sterile mixing decrease back in time: 

^i"»s= „0(7+,)3 ='"°»s(' + ^)''- (60) 
Here and 9g are parameters at the present epoch. The combination we plot in Figs. 1-3 

changes as 

sm^ggmg^ J J . 61 

(1 + zy 

The lines of constant induced mass sfiift witfi z to tfie left - to smaller masses and mixings. 
Tfiis means that in the past all cosmological bounds where satisfied. 

Already at the recombination epoch the mass of sterile neutrino becomes of the order 
10^ GeV for the present mass mPg ~ 1 eV and the mechanism of oscillation production does 
not work. Essentially, in this scenario the sterile neutrinos are not produced in the Early 
Universe and their concentration is negligible. The astrophysical and cosmological bounds 
we have discussed are not applicable. 

3) If S interacts with a massless or low-mass Majoron [86] 0, it can decay invisibly as 
5" — > or annihilate (see similar mechanism for the active neutrinos in [87]). If this decay 
is fast enough, T5 << 1 s, in principle, all astrophysical and cosmological bounds could be 
evaded. See [88] for the recent similar analysis. In the mass range ~ 1 keV such a fast 
decay can be achieved for the scalar coupling g ~ 10~^. 

In the low mass range some restrictions on the off-diagonal couplings gVaS^) can be 
obtained from their effect on the free-streeming condition for active neutrinos. The latter 
should manifest itself in the precision measurements of the CMB acoustic peaks. One can 
use these bounds also for sterile neutrinos. In [89] a limit on the active neutrino coupling 
gf < 1 • 10~^^(50 meV/ms)^ was obtained. According to Ref. [90], the limit is less severe and 
is absent if couplings with different active neutrinos are different, e.g., S has large coupling 
with Vt only. 



27 



5 Conclusions 



The main conclusions of this paper can be summarized as follows. 

1) Mixing of the active neutrinos with sterile neutrinos, singlets of the SM symmetry 
group, generates an induced mass matrix of active neutrinos which can be the origin of 
peculiar properties of the lepton mixing and neutrino mass spectrum. It opens an alternative 
possibility to understand possible new symmetries in the neutrino sector. In this way, one 
can explain the substantial difference of mixing patterns of quarks and leptons. 

Depending on masses and mixings of S, the induced active neutrino masses can be 
the origin of the dominant or sub-dominant structures of the neutrino mass matrix. For 
instance, the tri-bimaximal mixing can originate from the induced contribution. 

2) Apart from modification of the mass matrix of active neutrinos, there are direct mixing 
effects of S which can be observed in cosmology, astrophysics and laboratory experiments. 

The importance of the direct and induced effects depends on the range of S parameters 
considered. For ms > 300 MeV the induced effect dominates. The induced masses can 
reproduce the dominant structures of the active mass matrix. The direct mixing effects are 
neghgible. In the interval ms ~ (10~^ — 10^) keV the direct mixing effects dominate: the 
astrophysical and cosmological consequences of mixing are more important putting strong 
upper bounds on the induced mass. So, the latter can be neglected in the mass matrix of 
active neutrinos. In the narrow window ms ~ (0.1 — 0.3) eV the two effects are comparable. 
For the BBN bound ANi, < 1 the induced matrix can produce only small effects. If one 
additional neutrino is allowed by BBN, the induced masses can generate the sub-leading 
structures or even be comparable to the values of dominant mass matrix elements. 

3) New physics effects can relax or even lift the cosmological and astrophysical bounds 
thus opening a possibility to generate large induced matrix in the whole range of masses 
ms > 0.01 eV. Here interesting possibilities to notice are: cosmological scenarios with low 
reheating temperature, fast decay of S into a Majoron and neutrino, and the possibility of 
a soft mass ms which varies with time, a la Mass Varying Neutrino scenario. 
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Figure 1: The benchmark hnes Eqs. (33)- (35) and (55) versus the current astrophysical, 
cosmological and laboratory bounds on us — Ve mixing as described in the text. The colored 
regions are excluded in each case. The "thermalization" line and the two decay lines ts = Tree 
and Ts = Tjj are also shown. We show also the allowed regions for the warm dark matter 
and the LSND (3+1) as well as the region that could explain pulsar velocities. 
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Figure 2: Same as Fig.l but for — Vu mixing. 
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Figure 3: Same as Fig.l but for us ~ mixing. 
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Case 


Best Fit 

5 = \3 =0 


Experimentally allowed at 1 ct 

d = As = 


Free CP phases 


Normal 




' 3.2 6.0 0.6 ^ 

24.8 21.4 
V 30.7 ) 






' 2.5 - 5.0 2.7- 9.8 0. - 5.1 ^ 

19.9-30.3 18.1 - 22.9 
^ 24.5 - 34.0 y 






' 0.3- 5.0 0. - 10.8 0. - 11.1 ^ 

12.7 - 30.9 18.5 - 29.4 
^ 16.7-34.5 y 




Inverted 




' 48.0 2.8 3.7 ^ 

27.4 24.0 
V 21.7 y 






' 43.2 - 51.0 0. - 8.6 0. - 9.2 ' 

21.3-31.9 21.3- 25.6 
^ 17.8 - 28.2 ^ 


\ 
/ 




' 11.4-51.0 0.-39.0 0.-36.7 ^ 

0. - 32.1 4.6 - 26.7 
^ 0. - 28.2 ^ 




Degenerate 




' 200.0 0.5 0.4 ^ 

202.7 2.9 
^ 203.5 y 




/ 200.1 - 200.3 0.06 - 1.0 0. - 1.0 ^ 

202.1 - 203.6 2.4- 3.3 
\^ 202.5 - 204.1 y 






' 60.0 - 200.3 0. - 176.6 0. - 170.3 \ 

0.02 - 203.6 0.5 - 200.3 
^ 0.02-204.1 j 



Table 1: The reconstructed matrices for the normal and inverted mass hierarchies as well as the degenerate mass spectrum 
for mo = 0.2 eV. The moduli of elements are given in units of 1 meV. For each case we present: the matrix for the best fit 
values of the mixing parameters, the intervals of the matrix elements obtained varying the experimental values within 1 a 
for (5 = A3 = and also for free CP-violating phases. 
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